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Abstract. The Helfrich's shape equation of axisymmetric vesciles is studied. A sufficient condition on the 
physical parameters and some geometric properties are discovered for the formation of biconcave shape vesicles. 



The study of the biconcave-discoid shape of a 
red blood cell has been a continued interest in the 
last two decades. Theoretically, the shape of the cell 
is to minimize a certain bending energy. Like mini- 
mal surfaces and surfaces with constant mean cur- 
vature in geometry, the expected biconcave shape 
surface is a critical one for the bending energy func- 
tional. The early proposal for the bending energy 
by Canham [1] is purely geometric, which mainly 
involves the Willmore functional, [15, ch. 7]. It 
is well known among differential geometers that 
the unique minimum of the Willmore functional for 
topologically spherical vesicles is the round sphere, 
which differs from blood cells observed experimen- 
tally. Thus, Willmore functional or Canham's idea 
is not a good model for the shape of red blood cells. 
This is also observed by physicists [2]. 

It is known that the shape of blood cells and 
other biological membranes is closely related to 
the formation of lipid bilayer vesicle in aqueous 
medium. Based on this and the elasticity of lipid bi- 
layer, Helfrich proposed a modified bending energy, 
[4]. For a closed surface S, this bending energy is 
a combination of geometric quantities and physical 
parameters. The geometric quantities are the vol- 
ume V(S) enclosed by E, the area A(I]), the mean 
curvature H and the Gaussian curvature of S. 
According to the Gauss-Bonnet Theorem, the inte- 
gral of ii' is a topological constant. Within a certain 
topological class of E, Helfrich's bending energy can 
be reduced to, which we call Helfrich functional, 

= <j) {2H + cof dA + AA(S]) +pV(E) 

where cq, A — 2 A/ fee, andp — 2p/kc are the physical 
constant parameters of the functional interpreted 
as! IS the bending rigidity, cq the spontaneous 
curvature, A the tensile stress, and p — Po — Pi the 
osmotic pressure difference between the outer (po) 



and inner (pi) media. Here, we have taken a sign 
convention such that H is negative for spheres and 
Helfrich's numercial simulation produces a bicon- 
cave shape when cq is positive. 

The general minimizer of is not easy to find. 
In the past, much effort has been spent on study- 
ing axisymmetric solution to its variational equa- 
tion. This equation is usually referred to as shape 
equation of axisymmetric vesicles. There are nu- 
merous works on axisymmetric stationary vesicles, 
[2, 7, 14, 8, 12, 13, 9, 10]. It is stiU unknown to 
scientists what conditions on the physical parame- 
ters will guarantee a solution to the shape equation 
corresponding to a biconcave surface. 

In this letter, we provide a sufficent condition 
under which the Helfrich shape equation of axisym- 
metric vesicles has solutions of biconcave shape. Be- 
sides, we also report two necessary geometric con- 
ditions for the lipid bilayer vesicles, one for those 
with a reflection symmetry and another for those of 
biconcave shape. These conditions on the parame- 
ters are very mild but easy to verify. In particular, 
the case that cq > 0, A > 0, and p > is sufficient 
to ensure the existence of biconcave solution. The 
discussion on how the combination of the parame- 
ters Co, A, and p affects the existence of biconcave 
solution will be presented in another review paper 
in detail. 

The sufficient condition for the formation of bi- 
concave vesicle may be formulated in terms of a 
cubic polynomial, 

Q{t) ^t^ + 2cot^ + (eg + A)t - |. 

We proved that if all roots of Q{t) are posi- 
tive, then one can always find axisymmetric 
biconcave vesicles which are the stationary 
surface for the Helfrich functional. A pictorial 
summary is given at the end of this letter (p. 4). 
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For the necessary conditions, we observe that 
if Co > 0, the central curvature of an axisymmet- 
ric stationary vesicle of biconcave shape must be 
smaller than the first positive root of Q{t), more 
precisely, see (2). Furthermore, with the reflection 
symmetry orthogonal to the rotation symmetry, the 
Gaussian curvature along the "equator" (the inter- 
section circle of the axisymmetric vesicle with the 
plane of reflection) is explicitly given in terms of 
the radius rao of the "equator" and the polynomial 
Q{t), see (3). 

Equation and Requirements 

For an axisymmetric surface S, we take the rota- 
tional axis to be 2;-axis and the plane of reflection 
a;y-plane, then a biconcave shape is obtained by re- 
volving and reflecting a curve, given by a function 
z{r) > defined for r in some interval [0, Tqc], where 
Too is the radius of the "equator" . For a biconcave 
axisymmetric surface, we expect z{r) and its deriva- 
tive w = z' to have the following graphs. 




Usually, the Helfrich shape equation of axisym- 
metrc vesicles is written in term of the angle ip be- 
tween the surface tangent and the plane perpendic- 
ular to rotational axis [2, 6, 16]. For our discussion, 
it is convenient to rewrite it in terms of the deriva- 
tive of the graph w{r) — z'{r) — i&m}). Then the 
Helfrich shape equation of axisymmetrc vesicles be- 
comes an equation for w, which is 

/2 
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(1) 



The biconcave shape and symmetries are then 
translated into suitable initial and boundary con- 
ditions, such as, w{0) = 0; w'{0) = w'^ > 0; 
w{r) — oo as r — > r^; an integral condition that 

—007^ / w dr < 0; and w is required to have a 

unique local maximum and no other critical points. 



The lower order terms can be given by polyno- 
mials Q{t) or its quadratic part R{t) = Q{t) — t^. 
After multiplying with rw' , the equation (1) can be 
written as follow, 
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where K(r) 



w 



+ r^w'Q{K{r)). 
is the principal curvature 



of E in the meridinal direction. Indeed, 
K(r) — > w'q, as r — > 0, 



so the initial choice w'q is the central meridinal cur- 
vature or w'q^ the Gaussian curvature. Note that if 

Co = A = p = 0, equation (la) can be integrated 
directly to a solution z = z{r) which is a circular 
arc with radius 1 /w'q. This is exactly the case when 
the meridinal curvature is a constant. 



Necessary Conditions 

The lower order terms in (1) capture the changes 
of the meridinal curvature. Since Q{0) = R(Q) = 
—p/2 < 0, if is less than the smallest positive 
root of Q{t), then both Q{w'q) < and R{w'q) < 0. 
Our analysis shows that these negativity conditions 
arc essential for the meridinal curvature K(r) to 
change sign. Furthermore, if the physical parame- 
ters constitute a polynomial Q with only positive 
roots, a solution w{r) with w'q small enough al- 
ways satisfies all of our requirements. The threshold 
value of w'q can be calculated numerically if exact 
values of the parameters arc given. 

The meridinal curvature K{r) of the vesicle E 
is the most basic quantity in our analysis. From 
equation (lb), if cq > and the initial curvature at 
the center w'q is large, in the sense that R(w'q) > 0, 
then K{r) increases and a biconcave shape cannot 
be formed. So a necessary condition for formation 
of axisymmetric solution of biconcave shape in the 
case that cq > is 



R{w'q) = 2cot«^,' + {4 + ~X)w'q - I < 0. (2) 
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The following pictures show the behaviors of 

K(r) for two different values of with co = 1, 
A = 0.25, and p = I. For these parameter values, 
i?(0.277124) = 0. 




R{w'q) > 



R(w'f^) < 

The one on the left is taken with Wq = 0.278 so the 
necessary condition is not satisfied. In this case, k 
blows-up at finite radial distance and the solution 
only leads to a convex surface. The picture on the 
right is taken with w'q = 0.276. In this case, the 
necessary condition is satisfied, so k decreases and 
changes sign. But it only decays gradually and the 
integral area condition is not satisfied. As a result, 
the surface and its reflection intersects in M^. 

There is a more subtle boundary condition im- 
posed by the reflection symmetry. One needs to ver- 
ify that z{r) and —z{r) together form a stationary 
solution with respect to any variation near the plane 
of reflection. For the situation that we discuss, it 
turns out that this is true. Moreover, we obtain 
a relation between the radius Voo of the "equator" 
and the Gaussian curvature K{r oo) of the vesicle at 
any point on the "equator" , 




This relation is easily verified for the case that co = 
A = p = 0, in which, Q{t) = and K{roo)'^ = —. 

This is compatible with the fact that the solution 
surface is a round sphere. 

Sufficient condition 

To have a solution meeting all the requirements, 
we need to impose a stronger condition that Wq to 
be small enough, at least smaller than the positive 
roots of Q. This smallness assumption on w'q im- 
plies that not only R{w'q) < but also Q{wq) < 0. 
Both K(r) and w(r) start initially positive and be- 
come negative when r reaches a certain value ro 



(depending on the value of w'q, actually, compara- 

1/2 

ble to w'q ' when w'q is small). In fact, we discover 
that, for any point of infiection Tc of the graph of 
z{r) with K{rc) > 0, 

16 16 
— < lim < lim < — . 

6p w'g^O Wq w'g^O Wq P 

As a consequence of these inequalities, the graph of 
z{r) has a unique point of inflection. After r reaches 
To, with the assumption that all roots of Q{t) are 
positive, the lower order terms of equation (lb) are 
dominated and the function w{r) blows down mono- 
tonically to — oo in a finite radial distance. By fur- 
ther estimating the area according to w'q, the inte- 
gral area condition can be verified whenever w'q is 
sufficiently small. All these behaviors of the func- 
tions z(r) and w(r) with respect to Q and the choice 
of w'q are illustrated in our summary on p. 4. The 
mathematical details will be presented in a forth- 
coming paper. 

Discussions 

Finally, we would like to remark on other solutions 
which do not correspond to biconcave surfaces. 




The one on the left is a solution for cq = 1, A = — 3, 
p = 2 with w'q = 2. Note that w'q is greater than the 
positive roots of R{t) and Q{t). Moreover, Q{t) has 
negative roots and thus this is beyond the situation 
discussed above. In this case, after w = z' becomes 
negative, it may not decrease monotonically to neg- 
ative infinity, which is why a second "petal" may 
form. 

The other type is a rotational solution which is 
not symmetric under reflection in a:;y-plane [2, 3]. 
The picture is produced with co = — 1,A = — 1 
and p = \. The "upper" curve satisfles our equa- 
tion (1) but not the area condition. The "lower" 
curve and the reflected "upper" curve are both so- 
lutions locally, so there is a bifurcation of solution 
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at Too • The surface only closes up by an asymmetric 
one. It should be remarked that, from our numer- 
ical study, not all solutions in the upper quadrant 
may be closed up by an asymmetric counterpart. 
For example, if the upper part satisfies the area 
condition, the only possible lower part is the re- 
flection symmetric one. This suggests a possible 
uniqueness of biconcave solution. However, at the 
moment, there is not yet a mathematical proof. 
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Graphs of solution a 

(a) vj'q = 0.4 (b) ^u'^J = 0.25 


(c) w'g = 0.18 



Graph of Q and positions of w'q 
(co = 1, A = -0.3,p= 1) 




Graphs of solution z 

(a) »^ = 0.4 (b)TOj|=0.27 (c)m^=0.15 



Graph of Q and positions of w'q 
(co = -l, A = -0.25,p=l) 




Graphs of solution z 

(a) ^0 = 2 (b) ra' = 0.5 (c) vj' = 0.275 
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